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ABSTRACT 

When a planet falls onto the surface of its host star, the added high-metallicity 
material does not remain in the surface layers, as often assumed, but is diluted 
into the interior through fingering (thermohaline) convection. Until now, how- 
ever, the timescale over which this process happens remained very poorly con- 
strained. Using recently-measured turbulent mixing rates for fingering convec- 
tion, I provide reliable numerical and semi-analytical estimates for the rate at 
which the added heavy elements drain into the interior. I find that the rela- 
tive metallicity enhancement post-infall drops by a factor of ten over a timescale 
which depends only on the structure of the host star, and decreases very rapidly 
with increasing stellar mass (from about lGyr for a 1.3M Q star to lOMyr for a 
1.5M star). This result offers an elegant explanation to the lack of observed 
trend between metallicity and convection zone mass in planet-bearing stars. More 
crucially, it strongly suggests that the statistically-higher metallicity of planet- 
bearing stars must be of primordial origin. Finally, the fingering region is found to 
extend deeply into the star, a result which would provide a simple theoretical ex- 
planation to the measurements of higher lithium depletion rates in planet-bearing 
stars. 

Subject headings: hydrodynamics — instabilities — turbulence — planets and 
satellites: formation — stars: abundances 



1. Introduction 



Transit observations have recently revealed the existen ce of a number of very short period 
exoplanets ([Torres et al. I l2008h . Theoretical investigations (j Jackson et al.|[2009l : iLevrard et al. 



-2- 



20091 ) show that most of them are tidally unstable and doomed to fall onto the central star. 
These recent findings confirm the notion that, through disk migration fo llowed by tidal in- 
terac tions, the vast majority of planets formed may not actually survive (jLin fc Papaloizou 
19861 ). Back in the 1990s, this idea led Doug Lin to coin the colorful phrase "the last of the 
Mohicans" , to describe the few which - by one process or another - escape the fate of their 
siblings. 

The dire ctly-related question of whether the statistically-higher metallicity of planet- 
bearing stars ([Fischer fc Valentill2005l ) results from planetary infall or is of primordial origin 
remains unanswered to date. It is of crucial importance, however, since it bears on the debate 
between the core-accretion and gravitational-instability scenarios of planet formation. The 
various ideas proposed to find evidence in favor of the planetary infall theory are typically 
based on the premise that the added material is mixed with the outer convection zone 
only. The most commonly-discussed one relies on the very rapid decrease of the mass of the 
outer convection zone with stellar mass between 1M and 2M Q . The relative metallicity 
enhancement caused by the infall of similar- sized planets should there fore be much larger for 



has never been detected (ISantos et al.ll2001 



higher-mass stars than f or solar-mass stars ( Laughlin fc Adamslll997r) . This effect, however 



Fischer fc Valentill2005l ). a fact which may seem 



rather odd now that we know that infall events must regularly happen. Hence the question: 
what happened to the "other" Mohicans? 



Vauclairl ( 120041 ) was the first to question the original assumption that the added heavy- 



elements remain in the convection zone. She noted that, in the presence of adverse [i- 
gradients, a stably-stratified region could be destabilized and become turbulent through a 
well-known double-diffusive fingerin g instability , often called, by analogy with the oceanic 
context, "thermohaline convection" (ISternlll96Cf ). She argued that the induced mixing could 
slowly dilute the excess metallicity into the underlying radiative zone. 

Whether a system is stable or unstable to the fingerin g instability depends on the so- 
called density ratio R, defined in astrophysics (I Ulrichl 1 1 9 721 ) as: 



R 



V- V 



ad 



where, as usual, V — V 



and temperature, and /i is the mean molecular weight. The fingering instability takes place 
when 
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p and T are pressure 



1 < R < r 1 where r 



(2) 



kt is the thermal diffusivity and is the compositional diffusivity. A region with R < 1 is 
Ledoux-unstable to direct overturning convection. The fingering instability occurs in Ledoux- 
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stable systems, and is stabilized only when R > 1/r. Since 1/r is typically of the order of 
10 6 to 10 7 in stars, (1 Ulrichl 1 1 9 72l ; ISchmittl 119831). ey en a very modest adverse compositional 
gradient can trigger the instability. IVauclairl ( 120041 ) tentatively estimated that, through the 
added mixing, the decay timescale of the surface metallicity enhancement post-infall should 
be about 1,000 yrs, which would rule out the possibility of detection. At the time, however, 
fingering convection in the astrophysical context was much more poorly understood, and 
little was known about the strength of the turbulence actually induced. 



This, however, is no longer the case: using 3D numerical simulations, iTraxler et al. 



(120101 ) recently quantified turbulent transport by fingering convection in astrophysics. Their 
mixing laws are summarized in §21 I apply them, in £j3l to study the evolution of the 
/i— profile of a 1.4M star, polluted by a high-metallicity impactor. The results confirm most 
aspects of Vauclair's original idea but are now quantitatively reliable. In §H I propose very 
simple estimates of the process, which can be applied to a wide range of stars and impactors 
without the need for heavy-duty numerical integration. I conclude in £J5] by discussing the 
fundamental implications of these findings. 



2. Turbulent mixing by fingering convection 



Traxler et al.l ( 120101 ) found that compositional transport induced by fingering convection 

(3) 



is accurately modeled through the following diffusion coefficient: 
where v is the viscosity, 



R-l 
1/r- 1 



(4) 



is a reduced density ratic0, and g(r) is a universal function fitted from the numerical data: 

) ' if r e (0, 1) and g(r) =0 otherwise. (5) 



9{r) 



ioi e - 3 - 6r fi 



It is crucial to note that fl3]) has been obtained through experimental measurements, with 
no additional modeling, and contains no free parameters: it is a universal law for mixing 
by fingering convection. ITraxler et al.l ( 120101 ) also found that the induced tu rbulent heat 
trans port is negligible compared with the molecular and radiative contributions ( jDenissenkov 



20101 ). Hence it does not affect the background entropy profile of a star. 



1 R € [1,1/t] corresponds to r e [0,1]. The lower limit is now r = (below which the system is 
Ledoux- unstable) and upper limit is r — 1 (above which the system is fingering-stable). 
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In what follows, I adopt the following algorithm to calculate the total compositional 
diffusivity D^x.t): 

• For a given fi(x) profile (where x is the distance to the stellar center normalized by the 
stellar radius R+), calculate r. 

• If r > 1, set r to be 1 while if r < set r to be 0. 

• Set Dpfat) = D c ° nv (x,t) + D^ eT (x,t) + k^(x). 

In this last expression D c ° nv (x,t) is the fully convective diffusivity, chosen to be: 



rjconv 

D™ nv (x) = 



1 + tanh 



x — xi,(i) 
A 



(6) 



where x^ marks the transition between the Ledoux-unstable and the Ledoux-stable regions 
(i.e. below which r is strictly positive), A is the assumed thickness of the overshoot region 
(here = 0.001), and D conv is a very large convective value derived from mixing-length 

theory. 



3. Example of evolution of a polluted system 
3.1. The Model 

I now apply this new mixing law to model the evolution of the //-profile in a 1.4M Q star 
recently impacted by a y lMj planet. The backgroun d stellar model used was provided 



by P. Bodenheimer, see iGaraud &: Bodenheimerl (120101 ) for detail. The background mean 



molecular weight§| assumed is /i* ~ 0.61. Since the total mass of the outer convection zone 
of this star is significantly smaller than lMj, the metallicity of the polluted envelope is 
determined by that of the infalling planet. Here I choose a high- metallicity case with an 
initial molecular weight /io — 0.64 in the convection zone. 

In what follows, I compare the results of two related simulations, one in which fingering 
convection is ignored (by setting Z)^ nger to 0), and one in which it is taken into account. In 
both cases, all parameters, initial conditions and backgrounds are otherwise identical. The 
value of D conv (see (jUJ)) is taken to be equal to 10 8 cm 2 /s which is a few orders of magni- 
tude smaller than commonly-used convective diffusivities. This choice merely facilitates the 



2 I assume for simplicity that, in all objects, the Helium to Hydrogen abundance ratio is constant with 
Y/X = 0.39, and that the material is fully ionized, so that /i" 1 ~ (1.65 — 1.15Z). 
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numerical computation, as the timesteps would otherwise have to be prohibitively small. 
Indeed, it can be shown that the only effect of selecting a smaller convective diffusivity is to 
increase the timescale for the transient dynamical readjustment phase (see §1]). 



The initial /x-profile is: 



p(x) = //* + 



1 + tanh 



A 



(7) 



where x cz marks the base of the convection zone prior to infall. It then evolves according to: 



1 



d 



dp 

dt px 2 R\ dx 



where p is the background density, D^(x) is given by (EJ), where r is given by (jlj and R is 
rewritten more explicitly function of p(x) via: 



R(x) 



p{x)N 2 (x)R ir 



(9) 



where g(x) is the local gravity and A^ 2 = (g / H P )(V — V) with H p the pressure scaleheight. 
The functions p(x), g(x) and N(x), as well as v(x), k(x) and k^x) are extracted from the 
stellar model. 



Both simulations are evolved from the time of infall t = up to t = 700Myr. 



3.2. Results 

The evolu tion of the /^-profile goes through two phases, as originally discussed by 



Vauclairl (120041 ): a first rapid dynamical readjustment phase, and a second much slower 



fingering phase. These are now described in more detail. 

In the rapid adjustment phase (here t < 10,000yr), the depth of the outer convection 
zone first increases as the system becomes Ledoux-unstable. Convective mixing with the 
underlying material dilutes the added metals, and the surface mean molecular weight de- 
creases. As mixing proceeds, the compositional gradient at the base of the convection zone 
softens, and the convective zone retreats back to its original position, leaving behind a region 
of enhanced metallicity. After the end of this dynamical phase, the mean molecular weight 
of the convection zone is again more or less constant, but with a value p\ smaller than the 
original po. The region near the base of the convection zone is now marginally stable to 
the Ledoux-criterion. These points are illustrated in Figure [2j Note how up to now the 
difference between the models with and without fingering is practically indistinguishable. 
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Fig. 1. — Top panel: Evolution of fi(x) (equivalently Z/Z*) in the early dynamical adjust- 
ment phase of the outer convection zone. Shown are representative snapshots at (from top 
to bottom line in the convection zone) t — 0, t — 40yr, t = 160yr, t = 320yr, t = 640yr, 
1, 200, t = 2, 400yr, t = 4, 800yr and 10, OOOyr. The results of the simulations with fingering 
are indistinguishable from the one without. Bottom panel: early evolution of the position 
of the base of the convectively unstable region, x^. Its minimum is called x m j x in §4.11 The 
dashed line marks x cz of the unpolluted star. 



-7- 



Mixing by fingering convection begins to affect the system after about 10 5 years. The 
flux of high-/i material destabilizes the top of the stably-stratified region. The fingering 
region, which can be seen in Figure [3] as the extended interval where ~ 10 4 cm 2 /s, slowly 
but steadily deepens with timqf|. At these stellar parameters, the turbulent compositional 
transport is about a hundred times more efficient than molecular transport. By t ~ lOOMyr, 
the surface metallicity has nearly returned to its original value in the fingering case, but 
decays much more slowly in the absence of fingering. 



4. Estimated timescales 

The simulation presented above serves as a detailed illustration of the basic two-stage 
dilution process. However, one can also roughly approximate these two stages according to 
the following estimates, which can be easily applied to a wide class of stars without the need 
for numerical integration. 



4.1. Dynamical phase 

The end-product of the dynamical dilution phase can be modeled in the manner depicted 
in Figure HJ which is inspired from the results discussed in £j3j The convection zone has 
retreated back to its original position x cz , and has a mean molecular weight fi s = fi\. The 
compositional gradient in the region below has been smoothed out down to the lowest point 
of excursion of xl, which I call x mix . I assume that dfi/dx is constant between x mix and x cz , 
and therefore equal to (/ii — /^*)/(x cz — x m ix}- 

The convection zone dilution process conserves the total mass of added heavy elements 
so that, approximately, 

(fj.0 - /^)M CZ = (/ii - n*) mlx ^ cz , (10) 

where M cz is original mass of the outer convection zone, and M m j x is the mass of the spherical 
shell above x mix . This yields \i\ as: 

2M CZ 

-'"mix ~r i«cz 



3 Its progression would only be stopped where the background [i increases again in the nuclear burning 
region, not included in the computational domain x £ [0.5, 1]. 
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Fig. 2. — Evolution of the surface mean molecular weight fi s (t) of both simulations described 
in m 
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Fig. 3. — Evolution of fi(x) and D^(x) for the simulations described in £J3] Shown are 
times (from top to bottom line in the convection zone, top panel and from right to left in 
the fingering region, bottom panel) t = 10 5 yr, t = 10 6 yr, t = 10 7 yr, t = 10 8 yr and finally 
t = 7 x 10 8 yr. At t = 7 x 10 8 yr, the fingering region spans the entire computational domain. 
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The convection zone starts retreating when the system becomes Ledoux-stable, hence 
the marginal stability condition R = 1 uniquely determines x mix . At x mix , /i = /x*, and 
dfi/dx is roughly equal to a half the value it has in the region between x mix and x cz (since it 
is zero below x mix ). Using (|2J) and ffTTj) . the condition R = 1 becomes: 

/^(a;mix)A / ' 2 (^mix)-R* 2/z*iV 2 (a; mix ).R* 



fl f T .1^1 " f x . ) 2Mcz mo-M* 

9{X mix ) dx \ y^mix) M^+Mcic-iB 



12) 



This equation can be solved for x m ; x , which then yields /zi. For the 1.4M Q stellar model used 
here, the solution of f lT2|) is a; mix = 0.9165, so /ii = 0.636. This turns out to be a satisfactorily 
accurate estimate of the numerical results of §3] (where x mix = 0.919 and \i\ = 0.634, see 
Figure E]). 

Another test of this method was made integrating the evolution of the /i— profile of a 
1.3M star with /z* = 0.61, subject to a smaller initial enhancement (fio = 0.615) since its 
convection zone mass is now comparable with a Jupiter mass. The solution of equation (IT2"|) 
yields x m i x = 0.8695 and [L\ = 0.6147, which is again remarkably close to the numerical 
solution (x mix = 0.8705 and \l x = 0.6144). 



4.2. Fingering phase 

Given that even a very small adverse compositional gradient can destabilize the stratified 
background, the fingering region propagates rapidly into the interior. Moreover, within the 
bulk of this region, the turbulent compositional diffusivity remains close to its maximum 
value -D^tx = 101^/z/k^, as seen in Figure [3j Based on these observations, I model the 
dilution of the convection zone heavy-elements into the radiative interior in the following 
way. 

As in §4. H I consider a smoothly- mixed region, lying between x m ; x and x cz , where x m ; x is 
now assumed to deepen with time through a slow diffusion process. The evolution of x m i x (t) 
is derived by setting the diffusion timescale (via fingering) across the interval [x m i x , x cz ] equal 
to the time elapsed since infall, t: 



where a (which can be thought of as a spatial eigenvalue of this diffusion problem) is a 
number of order one which will be determined more precisely later, and D^^iL ls estimated 
from a typical value near x cz . As in §4.1[ by mass conservation: 

= + mJILio'" 1 " ' (14) 
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where /i s (£) is the evolving surface mean molecular weight, fii is its value at the end of the 
dynamical readjustment phase (estimated via the solution of ffTTT) and (fl2l . see §4.ip . and 
^mix(^) is the mass of the spherical region above x mix (t). 

The comparison of (fT4"l) with the exact numerical solution is shown in Figure [5j for the 
original run with the 1.4M star, and the additional test made with the 1.3M star (see 
§4.ip . The value of a is adjusted to fit the numerical solution for the 1.4M run: with 
a = 3.8 the approximate solution fits the numerical solution very well. The same value of a 
can then be used to estimate /x s (t) for the 1.3M run, and comparison with the true solution 
is again excellent. In both cases the main source of discrepancy turns out to be the estimate 
of fii, which can be off by about 10%. 

Having validated the semi-analytical estimates, I now use them on a variety of initial 
conditions and stellar masses (see Figure ED- The main result is that, when viewed in terms of 
the relative enhancement (/z s (£) — f-i*) I (/A) — the evolution of the surface mean molecular 
weight about lOMyr post-impact only depends on the structure of the host star and not on 
the initial conditions. The timescale for the relative enhancement to drop to 10% strongly 
decreases with increasing stellar mass, from about 600Myr for a 1.3M star, to 60Myr for 
a 1.4M star and 6Myr for a 1.5M star. Hence, one would be exceedingly lucky to detect 
any excess metallicity through planetary pollution on a 1.4M and 1.5M star. Meanwhile, 
higher-than-average-metallicity on a star of m ass 1.3M or s maller could be the sign of a 



recent (< lGyr old) event (e.g. HD 149026, see lLi et al.l (120081 )). but would of course require 



a very metal- rich, high-mass impactor to provide a detectably-high fi . 



Conclusion 



Thanks to the remarkable scientific contribution of iTraxler et al.l (120101 ). we now have 
reliable universal laws for turbulent mixing by fingering convection. In this Letter I apply 
them to the question of whether the metallicity enhancement caused by the recent infall of 
a planet onto its host star can be detected. 



As suggested by I Vauclairl (120041 ) . an infall event first causes a dynamical readjustment 
of the outer regions, triggered as the adverse compositional gradient destabilizes the system 
following the Ledoux criterion. This partially dilutes the added material into the upper 
radiative zone. Then follows a longer mixing phase, in which the radiative zone becomes 
unstable to fingering convection. The induced turbulence dilutes the remainder of the metal- 
licity excess into the radiative zone. The evolution of the surface metallicity in both phases 
can be modeled accurately through simple semi-analytical estimates provided in §4.11 and 
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I find that the relative metallicity enhancement decreases on a timescale which, after 
the first few Myrs, depends only on the stellar structure and is much shorter for higher-mass 
stars (~ lOMyr) than for lower-mass star s (~ lGyr), se e §4.21 for detail. This timescale 



is much longer than the one estimated by IVauclairl (120041 ). but is nevertheless quite short. 
This mass-dependence compensates the initial condition trend, whereby the post-impact 
metallicity is inversely related to the convection zone mass for a given impactor. These 
results explai n why no trend in metallicity vs. stellar m ass has been detected in planet- 
bearing stars (jSantos et al.ll200ll ; iFischer fc Valentil 120051 ). Furthermore, between the rapid 
dilution timescale for higher-mass stars and the the larger convection zone mass of lower- 
mass stars, the statistically-higher metallicity observed in planet-bearing stars cannot be due 
to planetary pollution: it has to be a primordial effect. 



Finally, note that by contrast with the assumption of I Vauclairl (120041 ). I find that the 
fingering region can extend very deeply into a star, and in particular, down to the Li- 
burning region. This would provide a si mple theoretical backing to the obse rvations of lower 
Li abundances in planet-bearing stars (lGonzalezll2008l ; llsraelian et all 120091 ). 
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Fig. 4. — Schematic of the simplified model. The dashed line shows the initial condition, and 
the solid line the current profile. The convection zone is homogeneous, with mean molecular 
weight n s . The radiative region adjacent to the convection zone is partially mixed, with a 
uniform gradient between x cz and x mix . In §4.11 /i, = ^i] in §4.21 /i c and x mix evolve with 
time. 
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Fig. 5. — Evolution of the relative enhancement in the surface mean molecular weight in the 
fingering phase. The symbols represents the numerical solutions, for the original 1.4M ex- 
ample described in $3]and for the 1.3M© test case described in §U The solid lines correspond 
to the estimated fi s {t) for each case and fit the data very well. Dashed lines correspond to 
runs with lower initial enhancements /j,q as shown. After about lOMyr, the evolution of the 
relative enhancement depends only on the stellar structure. Also shown as a dotted line is 
the estimated evolution of fi s in a 1.5M Q star. 



